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We show that modified gravity presents distinctive nonlinear features on the Cosmic Microwave Background 
(CMB) anisotropies comparing with General Relativity (GR). We calculate the contribution to the CMB non- 
Gaussianity from nonlinear Sachs-Wolfe effect in f(R) gravity and show that, contrary to GR's contribution 
which is typically < 0(1), the contribution in f(R) gravity is sensitive to the nonlinear structure of f(R) and 
can be large in principle. Optimistically, this gives an alternative origin for the possibly observed large CMB 
non-Gaussianities besides the primordial ones. On the other hand, such nonlinear features can be employed 
to provide a new cosmological test of f(R) or other modified theories of gravitation, which is unique and 
independent of previously known tests. 
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Understanding the origin and evolution of observed large- 
scale structure of the Universe plays a key role in cosmology. 
With the increasingly precise data of current and upcoming 
observations on cosmic perturbations (CMB, LSS), we are 
able to gain unprecedented understandings of both the Uni- 
verse and the fundamental physics relevant to it. An unique 
observable of cosmic perturbations is non-Gaussianity, which 
characterizes the nonlinearities of cosmic perturbations 13, Ht]. 
The possibly observed non-Gaussianities in CMB anisotropy 
have two main origins: the "primordial" non-Gaussianities 
of curvature perturbation £ yl Q], and the post-inflationary 
nonlinearities in the transfer from ( to anisotropies, or other 
secondary contributions 12J. The former is due to the inter- 
actions in the models or mechanism which generate the pri- 
mordial perturbations, and is thus expected to encode early 
Universe Physics. The later, however, which is due to the in- 
trinsic nonlinearities in gravitation and photon-baryon plasma, 
is often referred to as "contaminations" to the primordial non- 
Gaussianity, since the physics of CMB is well-understood and 
they are too small to be detected, especially compared with 
their possibly "large" primordial partners QUI]. 

In this Letter, we investigate the ability of these post- 
inflationary contributions to non-Gaussianities to probe gravi- 
tational theories alternative to General Relativity (GR). In the 
past several years, modified gravity theories (see e.g. n 1 Ofl for 
reviews) have gained much attention due to the unexpected 
later-time cosmic acceleration, which are employed as an al- 
ternative to dark energy models iflU [l2ll . Most of existed 
cosmological constraints to modified gravity theories (see e.g. 
1 19] for a recent review) come from the observation of CMB 
I HQ, LSS 0], weak lensing Q(J, solarsystem test Q3] 
as well as correction to the Newton's law IU8I1 . In this Let- 
ter, we propose a new observational window to test modified 
gravity, i.e. the nonlinear anisotropy due to gravitational per- 



turbations 10, 0, [U . As is well known, the large-scale radia- 
tion transfer comes from two aspects: the intrinsic anisotropy 
on the emission surface and the gravitational redshift under- 
gone by the photon from the emission surface to us. The later 
is kinetic and is irrelevant to the gravitational theory. How- 
ever, the former (as well as the nonlinear relation between 
metric perturbations and conserved curvature perturbation Q 
depends on the underlying gravitational theory. Thus, intu- 
itively, one may use this nonlinear transfer to probe the gravi- 
tational theories. In this Letter we focus on the simplest met- 
ric f(R) gravity. For the first time, we give the nonlinear 
generalization of the (ordinary) Sachs-Wolfe effect in f(R) 
gravity, up to the second order. Analogous to what happens in 
fc-inflation model with matter Lagrangian P(X, 4>) where pri- 
mordial non-Gaussianities are enhanced by the nonlinear de- 
pendence of P on X |@], in f(R) gravity the nonlinear trans- 
fer from primordial curvature perturbation into anisotropy is 
controlled by the nonlinear structure of f(R), and can also be 
enhanced in principle. 

We investigate a class of metric f(R) gravity with action 



S = 



(1) 



where we set 8irG = 1. The equation of motion for is 



I,rR^ - ;r/ (R) g M u + (g^o 



V M V„)/,j 



T, 
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where /,r = df /dR etc. and □ 

After decoupling, the CMB photon density remains as 
Planck distribution, which implies Tf/Ti = ujj/uji, where Tf 
(Tj) and o>/ (u>i) are the final (initial) temperature and photon 
energy respectively. The temperature anisotropy is thus given 
by 



AT 



T ~T 
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where lo q (u> e ) are scale factor and photon energy at the ob- 
server (emission surface). From (01 it is clear that the ob- 
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served anisotropy has two contributions: one is the gravi- 
tational redshift-induced anisotropy a u> /(a e ui e ), which is 
purely kinetic, the other is the intrinsic anisotropy T e /T e on 
the emission surface which depends on the dynamics during 
recombination and the gravitational theory. 
For later convenience, we introduce 



d\nH ^ _ d\nR _ d\nR' 
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where R' denotes dR/drj with 77 the comoving time. It is also 
useful to employ a set of dimensionless ratios, which charac- 
terize the nonlinear structure of f(R): 
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Recall that in Einstein gravity, we have e = 3(1 + w)/2 where 
w is the constant equation of state parameter, and e R = — 2e, 
e R < = 1 — 3e, (3 = 7 = <5 = 0. The "Compton parameter" B 
introduced in 111411 which was employed by several authors can 
be expressed as B — — /3e#/e, which should satisfies B > 
for the stability of linear perturbations ffl . 

We focus on the ordinary Sachs-Wolfe contributions, which 
dominate on the large scales. Thus we consider the scalar 
metric perturbations in Newtonian gauge 



ds 2 = a 2 (-e 2< W + e-™dx l dx l ) . 



(6) 



As mentioned before, the gravitational redshift contribution 
to the observed anisotropy is the same for all metric theory of 
gravitation, which is, up to the second order in $ J3>0-01 



-<J> + -<& 2 . 



(7) 



What is interesting is the intrinsic temperature anisotropy on 
the emission surface. Assuming the cosmic perturbation is 
adiabatic, in matter dominated era, the large-scale temper- 

ature fluctuation has the relation: T e /T e = (p 1 /p 1 ) i = 
_ 1 

(pm/Pm) 3 ■ To determine the matter density perturbation, we 
have to make use of the generalized Einstein equation @. In 
GR, on large scales it yields a non-perturbative relation |@]: 
Pm/pm = e -2 *. In f(R) gravity, up to the second-order in 
metric perturbations, we find 
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(2-a 1 )$+i(a 2 -a 2 )$ 2 , (8) 



where a\ and a 2 are combinations of derivatives of / with 
respect to R, 
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(9) 



1 + /j(e - 1) + (2/3 + j)e R 
2e - - pe R (e R , - 2) ' 

10 /3-l + 7(e-l) + (37 + *)efl , lm 
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where parameters are defined in dUi-©. It is useful to note that 
in GR, o"i = (72 = 0, thus the familiar relation \n(p m / p m ) = 



—2$ is recovered (up to the second order). The intrinsic tem- 
perature fluctuation on emission surface now takes the form: 



Te 



1 



1 



1 + - (-2 + tn) $ + — (4 - 4 CTl - 2a\ + 3a 2 ) $ 2 . 



(11) 

Finally, collecting (O and (fTTT i. after some manipulations, 
we get the nonlinear Sachs-Wolfe contribution to the large 
anisotropy up to the second-order in $ in f(R) gravity: 

AT 1 1 

— = - (1 + ax) (1 + 2a! - 2aj + 3<7 2 ) $ 2 . (12) 

In the above, quantities are evaluated at the emission surface. 
([T2I 1 generalizes the familiar results Of- = e*/ 3 in GR |@]. 

Our next goal is to relate the observed temperature 
anisotropy to the primordial curvature perturbation (, which 
encodes the information in the very early universe and is 
the most frequently used variable in evaluating the primor- 
dial non-Gaussianities in the literature. It is well known that 
on large scales and for adiabatic perturbation, there is a non- 
perturbative and conserved quantity which can be identified 
as nonlinear curvature perturbation in uniform-density slices 
ll2lll (see also 1I22I1 ). defined as 



T 1 r d p 

C = — * + - / — . 



(13) 



The conservation of ( relies on the fact that the energy- 
momentum tensor of the perfect fluid is conserved, and is ir- 
relevant to the gravitational theory. Our next task is to set the 
initial conditions for $ and on the emission surface up to 
the second-order in £. In matter-dominated era (p = 0), (fT3l l 
can be integrated to give 



3 Pm 



(14) 



Together with (|8), it will gives ( in terms of metric perturba- 
tions $ and *: ( = <£]. 

The next step is to determine the constraint between $ and 
in f(R) gravity. As in GR, this can be done by employing 
the spatial traceless part of the generalized Einstein equation. 
Applying didj on both sides of SjJ = T^ F , where Sy is 
defined in (O, we are able to solve ^ in terms of $ on large 
scales, which gives, at linear order 



*i = (1-2/3)$, 



(15) 



where f3 is defined in (0. ( fT3T > is the well-known result, which 
implies that in f(R) gravity, $ 7^ $ even at the linear order 
for perfect fluid. At the second-order, we find 

* 2 = A' 2 [$] = d- 4 [(3A-2^ + 8/3 2 +4 7 ) («9 2 $) 2 
+ (A + 6£ - 4/3 2 + 87) {d t d^f 
+4 (A + 2/3 + 4 7 )a t $a t 9 2 $ 
+4(/3-/3 2 + 7 )$a 4 $], (16) 

with 

2(l-/3 + (3/3/2 + 7) ej i) 2 
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2e + pe R (2 - e R ,) - 



(17) 
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where various parameters are defined in ©-(O. (fT5Tl-(fT6]l 
generalize the corresponding constraints in GR |2|, l8[ . Note 
the above parameters are all evaluated on the background. The 
GR results corresponding to (TT3T> - (TT6b can be easily recovered 
after plugging the GR values of parameters. Having these re- 
lations in hand, now it is straightforward to get the initial con- 
ditions for $ up to the second order in £: 



$1 



3( 



5 - 6/3 - a x ' 
9 [(q 2 - a\) C 2 - &K 2 [C]] 
2 (5 - 6/3 - aif 



(18) 
(19) 



where K2 [•] is defined as in ( TToT ). 

From the nonlinear Sachs-Wolfe expression (TTZb and the 
exact expression for curvature perturbation in matter era (TPfl) . 
together with (TT~8b-(fT9T>. a nonlinear mapping from the con- 
served curvature perturbation to the temperature anisotropy 
can be got: = + ( 4+1 H , which at the 



linear order is 



• AT\ 
. T 1(2) 



J /(I) 



l + (Tl 
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and at the second-order 
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where 



b{pi,P2) =b - hg (pi,p 2 ) 



(21) 



(22) 



other only in the different coefficients &o an d &i- In General 
Relativity, it is easy to verify that J2, Htl 



b(px,p 2 ) 



GR 
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(26) 



with the same momentum factor g(px,P2)- The GR expres- 
sion d26l i implies that, due to the smallness of the numerical 
coefficients, the nonlinear mapping (f2Tb would not contributes 
significantly, especially in the so-called "squeezed" momenta 
configurations of non-Gaussianity |@], which is a long well- 
know conclusion in GR. In f(R) gravity, however, the coeffi- 
cients bo and b\ depend on the nonlinear functional structure 
of f(R), and thus can be numerically large in principle. In 
Figfl] we plot the dependence of bo and b\ on parameters j3, 
7 for a ACDM expansion history for illustrative purpose. 
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with momentum-independent coefficients 
1 



b 



(5-6/3-(7i) J 



5-6^(7 + 2(4-cri) t r 1 +3ct 2 
+36 (/3 2 - 7 )(l + ai) 

+9(71 - O- 2 (15 + (7l) + 18f7 2 



bx 



3 (l + ai)(3A + 6/3(2/3- 1)) 
2 (5 - 6/3- cri ) 3 



(23) 
(24) 



and momentum-dependent factor 



9{P\,P2) = 1 + 2 



Pi 



(pi +P2) 



3^%. (25) 
(Pi +P2) 



Note in (I23b-(l24b. parameters are evaluated at the time of 
emission (e.g. decoupling). The above expression has been 
arranged in such form so that g(pi,P2) —> as px —> 
or P2 — > 0, which denotes the so-called "squeezed" limit 
of momenta configuration for the bispectrum. Thus in this 
"squeezed limit", b(px,k — pi) — > bo- A large bo implies 
a large contribution to the "local"-type non-Gaussianity. The 
functional form of &(pi,p 2 ) is rather general, the contribu- 
tions in Einstein gravity and f(R) gravity differ with each 



FIG. 1: (Color online). 60 and 61 as functions of j3 with diverse 
values of 7 and vanishing 5 for the ACDM expansion history with 
SIa = 0.74. In the matter-dominated era, the parameters introduced 
in © are e = 1.5, en = —3 and e R i = —3.5 respectively. We 
assume the range of values of /3, 7 and 5 ensures such an expansion 
history. 

Most efforts of investigating the non-Gaussianity focus on 
the primordial ones, since they are expected to encode the 
early Universe physics. In this Letter, we explore the abil- 
ity of "non-primordial" (more precisely, "post-inflationary") 
non-Gaussianity to probe gravitational theories. In the con- 
text of f(R) gravity, we calculate the nonlinear CMB tem- 
perature anisotropy in the Sachs-Wolfe realm. We derive the 
explicit nonlinear mapping from the primordial curvature per- 
turbation to the observed anisotropy, up to the second order. 
We find that the amplitude of this nonlinear mapping, which 
is typically < 0(1) in GR H H, is controlled by the 
functional structure of f(R) and thus can be large in prin- 
ciple, as the primordial non-Gaussianities can be enhanced by 
the non-canonic ality of kinetic term(s) of inflaton field(s). In 
this Letter we focus on the Sachs-Wolfe contributions which 
dominate on the large-scales. A full treatment of the nonlin- 
ear transfer on all scales involves solving the full Boltzmann 
equations l2(ill . which we leave for future investigations. 
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The result presented in this Letter has at least two impli- 
cations. On one hand, it gives an alternative origin to the 
possibly observed non-Gaussianity — especially the so-called 
"local type" — beyond primordial ones. On the other hand 
and might be the most interesting, this sensitivity of nonlin- 
ear transfer to the structure of f(R) can be employed to pro- 
vide a new test of f(R) gravity, which is unique and indepen- 
dent of previously known tests, e.g. from spectra of CMB 
or LSS lfl3l - [l6ll . As has been emphasized, the momenta- 
dependent factor (l25l l is universal for f(R) theories but is 
also distinctive comparing to momentum shape functions of 
primordial bispectrum ]6[]. Thus, by constructing appropri- 
ate templates, we are able to disentangle f(R) contributions 



from other primordial/non-primordial sources in the finally 
observed non-Gaussian signals. In light of current and upcom- 
ing observational constraints on non-Gaussianities, we may 
have already been able put new constraints on f(R) or more 
general modified theories of gravitation. 
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